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Abstract 

Motivated by the isotropy of the CMB spectrum, all existing studies of magnetised cos- 
mological perturbations employ FRW backgrounds. However, it is important, to know the 
limits of this approximation and the effects one loses by neglecting the anisotropy of the back- 
ground magnetic field. We develop a new treatment, which fully incorporates the anisotropic 
magnetic effects by allowing for a Bianchi I background universe. The anisotropy of the 
unperturbed model facilitates the closer study of the coupling between magnetism and ge- 
ometry The latter leads to a curvature stress, which accelerates positively curved perturbed 
regions and balances the effect of magnetic pressure gradients on matter condensations. We 
argue that the tension carried along the magnetic force-lines is the reason behind these 
magneto-curvature effects. For a relatively weak field, we also compare to the results of the 
almost-FRW approach. We find that some of the effects identified by the FRW treatment 
are in fact direction dependent, where the key direction is that of the background magnetic 
field vector. Nevertheless, the FRW-based approach to magnetised cosmological perturba- 
tions remains an accurate approximation, particularly on large scales, when one looks at 
the lowest order magnetic impact on gravitational collapse. On small scales however, the 
accuracy of the perturbed Friedmann framework may be compromised by extra shear effects. 

PACS numbers: 9880H, 0440N, 4775, 9530Q, 9862E, 0420 



1 Introduction 

Over the years, the implications of primordial magnetic fields for the formation and the evolution 
of the observed structure have been the subject of continuous theoretical investigation. The 
majority of theses studies are Newtonian, or semi-relativistic, and all of them are confined within 
perturbed Friedmann-Robertson- Walker (FRW) models. Mathematically speaking, however, 
the spatial isotropy of the Friedmann universe is not compatible with the presence of large-scale 
magnetic fields. Nevertheless, one may still preserve the FRW symmetries by postulating a 
random background field (e.g. [Q, or by allowing for a uniform field that is too weak to 
destroy the Robertson- Walker isotropy (e.g. |J-[||]). In either approach the anisotropy induced 
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by the magnetic field to the zero-order stress tensor of the universe is treated as a perturbation. 
Both treatments, apart from their different initial set up, proceed essentially in the same way 
and they should reach the same conclusions. 

Strong motivation for adopting an FRW background comes from observation. Current mea- 
surements of the Cosmic Microwave Background (CMB) anisotropy suggest an energy density 
for the primeval field much smaller than that of matter at recombination j6j. Similar limits, 
though less stringent, are placed at nucleosynthesis by Helium-4 measurements @. One may 
also use dissipative arguments to limit the strength of primordial magnetic fields on small scales 
H . Intuitively, the impact of such a weak field on the isotropy of the smooth FRW spacetime 
should be negligible. Nevertheless, pending on more detailed treatments, one might still ask: 

• What is the exact nature of the effects induced by the anisotropic magnetic stresses? 

• How important are the deviations from the results based on perturbed FRW models? 

We address these questions by employing a Bianchi I background, the simplest anisotropic 
cosmology that can naturally accommodate large-scale magnetic fields. We also use the covariant 
Lagrangian perturbation formalism Q, which has proved powerful in dealing with issues of 
similar complexity in the past. This article extends the work of |B|], where the exact equations 
of a general magnetised cosmology were derived and subsequently linearised about an FRW 
spacetime. Here, we linearise the same equations about a Bianchi I universe supplied with a 
large-scale magnetic field and a perfectly conducting ideal fluid. Our main aim is to identify 
and analyse the kinematic and dynamic effects of the anisotropic magnetic stresses. In the 
process we will also test the accuracy of the FRW-based treatments of magnetised cosmological 
perturbations. 

We provide the full set of the linear equations and then isolate the effects of the background 
magnetic stresses, by treating the field as an eigenvector of the shear tensor to zero order. This 
way we can identify the deviations from the almost-FRW results that are caused purely by 
the background magnetic stresses. Note that background anisotropics that are not related to 
the magnetic presence are not relevant to our study and therefore not included. At the linear 
level, however, the field is not assumed to be a shear eigenvector. Thus, in the perturbed 
universe the shear is no longer the pure result of the magnetic stresses. We consider a universe 
filled with a single perfectly conducting barotropic fluid and compare to the predictions of 
the FRW-based treatments. Generally speaking, the impact of a cosmological magnetic field 
has an isotropic as well as an anisotropic component. The former results from the magnetic 
contributions to the total energy density and isotropic pressure, whereas the latter is due to 
the magnetic stresses. Our analysis shows that although the Friedmannian limit identifies the 
isotropic component exactly, it can only address the 'average' anisotropic effects of the field. 
Note that the 'average' is taken over the three linearly independent spatial directions rather than 
over the whole space. Thus, as one should expect, the isotropy of the background spacetime sets 
certain limits to the accuracy of the standard treatment. In this report we show that some of 
the magnetic effects identified by the FRW approximation are in fact direction dependent, with 
the key direction determined by the zero-order field vector. Even when the field is relatively 
weak, certain quantities, such as the curvature, continue to contribute differently parallel to the 
magnetic force-lines than orthogonal to them. Clearly, in the linear regime, such directional 
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dependence results from the background anisotropy and it is not accessible through an FRW 
background. 

Of particular interest are the magnetohydrodymanical effects on the average volume ex- 
pansion of the universe. We find that density perturbations co-linear to the background field 
contribute more efficiently to the deceleration rate. Also, only spatial curvature perturbations 
parallel to the zero-order magnetic vector affect the volume expansion. The latter effect is of 
geometrical origin, resulting from the vectorial nature of the magnetic field. It was first iden- 
tified in [[!(]], although there the isotropy of the FRW background meant that the effect was 
independent of direction. Here, as in |uj, the magneto-curvature term in the Raychaudhuri 
equation decelerates negatively curved regions, but accelerates those with positive local spatial 
curvature. The crucial difference made by the background anisotropy, is that it reveals the 
prominent role played by the field lines themeselves. We can now attribute the unconventional 
magneto-geometrical effects to the negative pressure, the tension, experienced along the mag- 
netic lines of force. This property of the field, together with the coupling between magnetism 
and geometry, gives rice to a curvature stress. The latter reacts to the effect of local magnetic 
pressure gradients by modifying the expansion rate of the region accordingly. This relativistic 
effect is analogous to the classical curvature stress exerted by field lines with a nonzero local 
radius of curvature. 

Direction dependent contributions are also identified in the evolution equations of the density 
gradients. Curvature effects, for example, are again confined in the direction of the background 
field. At the weak field limit and on super-horizon scales, this change in the curvature contri- 
bution is what separates the Bianchi I perturbation equations from their FRW counterparts. 
Nevertheless, such directional effects can only alter the lowest order magnetic impact during 
the radiation era. Even then, the predictions of the FRW approximation still hold, at least in 
certain critical cases. Our results verify the point made in [11], that magneto-curvature effects 
tend to counterbalance the 'pure' magnetic impact on the evolution of density gradients. In 
particular, for zero curvature input the field was found to suppress slightly the growth of the 
density contrast ||10||, in agreement with the Newtonian treatment of ||. However, when cur- 



vature contributions were included in the equations, this inhibiting effect vanished [11]. Here, 
we verify these statements and argue that the magneto-geometric contributions can also reverse 
the pure magnetic impact on gravitational collapse. When the curvature perturbation along 
the field lines is maximum, we find that the density contrast grows faster than in standard 
non-magnetised cosmological models. Thus, the combined magneto-curvature action opposes 
the effect of the pure field on local matter condensations. Again, the explanation comes from 
the tension experienced along force-lines of the field. When coupled to geometry, this magnetic 
property leads to a curvature stress, which always tries to balance the effect of the magnetic 
pressure gradients. 

On sub- horizon scales, we find that curvature complexities are accompanied by similar di- 
rection dependent shear effects, even when the field is weak. These shear contributions result 
from the kinematical implications of the background magnetic stresses, and do not appear in an 
almost-FRW treatment. One might be able to invoke fluid 'backreaction' arguments and ignore, 
to first order, the aforementioned small-scale shear effects. In this case, one recovers the stan- 
dard FRW linear equations, provided that the curvature perturbation is averaged. Nevertheless, 
one should be rather careful when discussing the small-scale magnetic effects within perturbed 
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FRW cosmologies. In contrast, super-horizon weakly magnetised density perturbations proceed 
basically as predicted by the Friedmann approximation. 

This paper is organised as follows: Section 2 provides a brief review of the covariant and 
gauge-invariant approach to cosmology with references for further reading. In section 3 we 
linearise the exact equations about a magnetised Bianchi I background. Particular properties of 
the unperturbed model are discussed in section 4. In section 5 we apply the linear equations to 
the case of a barotropic medium, and examine the field impact on the model's kinematic and 
dynamic evolution. Section 6 addresses the weak field limit and section 7 discusses solutions of 
the perturbation equations. Finally, we summarise our conclusions in section 8. 



2 The Covariant Formalism 



In the covariant Lagrangian approach to cosmology (see \12\ for an up-to-date review) one makes 
a physical choice for the 4-velocity u a of the matter. Then, one uses h ab = g a b — u aUb, where g a b 
is the spacetime metric, to project into the instantaneous rest spaces of the comoving observers.^] 
Note that these rest spaces form an integrable 3-dimensional hypersurface only for zero vorticity. 
Nevertheless, we will keep referring to them as 'spatial' or '3-spaces' for simplicity. Thus, the 
fluid 4-velocity introduces a unique 1+3 splitting of the spacetime into time and space. We 
can now decompose tensor fields and tensor equations into their timelike and spacelike parts. 
Without going into much detail, we provide some key decompositions essential to our analysis. 



2.1 The Basic Equations 

The kinematical properties of the motion arise from splitting the covariant derivative of u a into 
its irreducible parts 

V b u a = \®h ab + o ah + e abc uj c - u a u b , (1) 

where e abc is the rest-space volume element, B = V a ii a is the rate of the volume expansion, 
a a b = V \a u b) ls t ne shear, u> a = — ^curru a is the vorticity and ii a = u b V b u a is the acceleration.^] 
The volume expansion defines an average length scale S via S/S = |©.0 

The medium is characterised by its energy-momentum tensor. Relative to a comoving ob- 
server, the stress tensor of an imperfect fluid decomposes as 

T ab = H*u a u b + p*h ab + 2q* a u b) + vr* b , (2) 
1 We employ geometrised units with c = 1 = 8nG and our metric signature is (— H — h +). 

2 Following |is||, we use angled brackets for the projected, symmetric, trace-free (PSTF) part of tensors and for 
the orthogonal projections of vectors. The same notation also represents orthogonally projected covariant time 
derivatives. 

3 For comparison reasons, the current article adopts the notation of jHJ with a number of changes motivated 
by jl2| and JL3| |. In particular, totally projected covariant derivatives are now represented by V instead of '^V, 
II a (, has replaced M a t as the anisotropic magnetic pressure tensor, the 3-Ricci tensor and the 3-Ricci scalar have 
changed from ^ R a b and K to 7L a b and 1Z respectively, while the scalar B represents now the dimensionless ratio 
S 2 X7 2 H 2 1 H 2 rather than the Laplacian X7 2 H 2 . Other changes include a difference in the definition of the Alfven 
speed (see §4.2.2) and the use of the vorticity vector ui a instead of the antisymmetric tensor uj a t- Finally, we 
employ the 'streamline' definitions for the covariant curls of spacelike vectors and tensors (see 113]]). 
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where fx* and p* are respectively the relativistic energy density and the isotropic pressure, q* 
is the energy flux vector and 7r* b is the PSTF tensor representing anisotropic pressures. The 
detailed physics is encoded in the equations of state. For a perfect fluid, with zero flux and 
anisotropic stresses, Eq. @ reduces to 

Tab = fJ>u a u b + ph ab . (3) 

Similarly, in the absence of electric fields, the stress tensor of a pure magnetic field H a splits 

as 

T ab = \H 2 u a u b + \H 2 h ab + U ab , (4) 

where H 2 = H a H a . The latter provides a measure of the magnetic energy density and isotropic 
pressure, while 

n a b = ^H 2 h ab — H a H b = —H( a H b } (5) 

is the PSTF tensor that represents the anisotropic magnetic stresses. 

The Weyl tensor C a bcd represents the long-range gravitational field, namely tidal forces and 
gravity waves. Together with the Ricci tensor R ab , which is decided locally by matter, it de- 
termines the Riemann curvature tensor R a bcd completely. For a comoving observer we have 


C ab cd = 4 (u [a u^ c + h { J c ) E b] d i + 2e abe u^H^ + 2u [a H b]e e cde , (6) 

where E ab = C acbl ju c u d and H ab = ^£ a cdC cd be u e are respectively the 'electric' and 'magnetic' 
parts of C abcd - 



2.2 The Basic Gauge- invariant Variables 

The covariant Lagrangian formalism is the foundation for a covariant and gauge-invariant pertur- 
bation theory |J, which provides an alternative to the metric-based gauge-invariant formalism 
|15|| . In the covariant approach, one deals with the inhomogeneous and anisotropic spacetime 
directly, instead of perturbing away from the smooth background universe. Thus, one obtains 
the exact (fully non-linear) equations first, before linearising them about a chosen background. 
This is clearly a major advantage, because it allows one to address cosmological models more 
general than the perturbed FRW universes. Indeed, the covariant approach has already been 
employed to study perturbed Bianchi I cosmologies filled with a non-magnetised perfect fluid 
fli~6l| . Here, we use the same formalism to analyse, for the first time, almost Bianchi I spacetimes 
permeated by a large-scale magnetic field. 



The covariant approach utilises the Stewart- Walker lemma [17j, to define a set of gauge- 
independent tensors. These describe inhomogeneities in the key quantities and have a trans- 
parent geometrical and physical interpretation. For our purposes, the basic gauge-invariant 
variables are || 

V a = -SV af ,, Z a = SV a Q and M ab = SV b H a . (7) 
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The above respectively represent spacelike variations in the energy density of the fluid, the 
volume expansion and the magnetic field vector, as these are seen by a pair of neighbour- 
ing fundamental observers^] They vanish in an exact Bianchi I spacetime, thus satisfying the 
Stewart- Walker criterion for gauge-invariance. Note that gradients in the fluid pressure are rep- 
resented by the gauge-independent variable Y a = V a p. For a barotropic fluid, however, Y a is 
always given in terms of T> a . 

3 Perturbed Magnetised Bianchi I Universes 

The non-linear equations that describe density perturbations of a perfectly conducting medium 
in the presence of a large-scale magnetic field have been derived in ]5|. Here, we linearise these 
equations about a Bianchi I spacetime. Covariantly, the dynamics of the model is characterised 
by the energy density and the pressure (p) of the matter. The magnetic presence also brings 
into play the field's energy density and isotropic pressure (both expressed via H 2 ), as well as 
anisotropic stresses represented by H a b- Kinematically, the Bianchi I cosmology is determined by 
the average expansion scalar (B) and by the shear tensor (a ab ), which describes departures from 
the isotropic expansion. Anisotropies in the spacetime curvature propagate via the electric Weyl 
tensor (E ab ), while the spatial sections of the Bianchi I spacetime are flat. During linearisation 
all these variables are treated as zero-order. First-order quantities are those that vanish in 
the background and therefore satisfy the requirement for gauge invariance [jT^| . These are the 
acceleration (u a ), the vorticity (u; a ), all variables representing spatial curvature, the magnetic 
Weyl tensor (H ab ) and the inhomogeneity variables V a , Y a , Z a and M. ab . Terms of higher 
perturbative order will be disregarded. 

3.1 The Medium 

According to Eqs. (||) and (Q), the energy-momentum tensor for a magnetised spacetime filled 
with a single perfectly conducting ideal fluid is 

T ab = (m + \H 2 ) u a u b + (p + \H 2 ) h ab + U ab . (8) 

where the assumption of infinite conductivity ensures the absence of electric fields || . Comparing 
Eqs. (||) and (|2|) we see that our magnetised medium corresponds to an imperfect fluid with 
/i* = fi + ^H 2 , p* = p+^H 2 , q a = and ir* b = H ab = —H^ a H b y Thus, magnetism contributes to 
the total energy density and isotropic pressure, while it is solely responsible for the anisotropic 
stresses. Note that zero electric field implies a vanishing Poynting vector, which in turn explains 
the absence of a energy flux vector in Eq. (|8|). 

3.2 The Linear Equations 

Next, we provide the key relations that describe the linear evolution of our magnetised Bianchi I 
cosmology. Unless stated otherwise, all equations derive from the non-linear formulae given in 

4 in a perturbed Bianchi I spacetime the aforementioned geometrical interpretation of the key inhomogeneity 
variables holds only for weak shear (see Appendix). 
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[]|].[] In particular, we employ a combination of propagation and constraint equations. These 
are: 

(i) The fluid conservation laws, represented by the continuity equation 

ii = -n(i + w)e (9) 

and by Euler's formula 

ti(l + w + $h)ua = -Y a - e abc H b cm\H c - u b U ba , (10) 
where w = p/fi, h = H 2 / u and cm\H a = e abc V b H c . Generally, w is allowed to vary with 

w = -(i+w)(cl-w)e, (ii) 

where c 2 = p/fi. Therefore, when w = 0, then c 2 = w = constant. Note the total absence 
of magnetic terms in Eq. (|9|), since the field's energy density is conserved separately as a 
consequence of Maxwell's equations. 

(ii) The kinematic propagation equations, comprising Raychaudhuri's formula 

9 = -\Q 2 -i/x(l + 3iy + h) -2a 2 + A + A, (12) 

where A = V a u a (to first order) and A is the cosmological constant, the vorticity propagation 
equation 

tu a = -§0o/ a + a ab io b - ^curlu a (13) 
and the shear propagation equation 

<j(ab) = -3@0"a6 - Vc{aV C b) + V( a ttfe) + \li a b ~ E ab . (14) 

The above are supplemented by an equal number of constraints. In the nomenclature of []Ej] , 
these are the (Oa)-equation 

|V a 6 = V b a ab + curla; a , (15) 

the vorticity-divergence constraint 

V a uj a = (16) 

and the -H^-equation 

H ab = (curia) ab + V( a a; b ) , (17) 

where (curls')^ = £ c d(aS c S d b) f° r every PSTF tensor S ab [13]. Note that constraint ( p|) ensures 
that the uj a is a solenoidal vector to first order. 

(iii) The magnetic equations, obtained from the covariant Maxwell equations [19]. For a 
perfectly conducting medium, the latter decompose into one propagation equation 

H {a) = -l&H a + a ab H b + e abc H b oj c , (18) 



5 For a quick cross-check, the reader may compare our equations to those of Sec. 2 in . The latter must be 
adapted to the particular imperfect fluid defined by Eq. (H) . 
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namely the magnetic induction equation, and three constraints 

e abc u b H c + curlfla = J {a) , (19) 

2uj a H a = e , (20) 

V a H a = 0, (21) 

where is the projected current density and e = —J a u a is the charge density. According to 
Eq. (^o|), an overall electrical neutrality (i.e. e = 0) ensures that the vorticity vector and the 
magnetic field are orthogonal. Also, Eq. ( |2l| ) guarantees that H a is a solenoidal. Contracted 
with H a , the induction equation provides the conservation law of the magnetic energy density 

(h 2 ) = -±QH 2 -2o ab Tl ab , (22) 

where the background anisotropy has added the last term of the right hand side. Thus, unlike 
perturbed FRW models, the magnetic energy density no longer drops as S~ 4 . Also, from Eq. 
( |27| ) bellow we obtain the expression 

n a fe = §6<7 a fc — 2<7 c ( a <7 c fe } — 2E ab + TZ^ a b) , (23) 

that connects Tl ab to the rest of the linear anisotropic sources. 

Note that, so far, all equations have been derived from the fully non-linear relations given 
in ||. Next we provide a new pair of linear equations that govern the behaviour of H a b- In 
particular, the time derivative of (|5|) and Eqs. ( |lS| ) and (22) lead to the propagation equation 
of the anisotropic magnetic stresses 

n<afc> = -f en afe + 2II c<a (a c b) + e cd b) uj d ) - \iiha ab . (24) 

Similarly, by projecting the gradient of @ orthogonal to the fluid flow we arrive at the constraint 

V b n ab = e abc H b C ui\H c - lV a H 2 . (25) 

(iv) The equations that determine the geometry of the observer's rest space. The latter is 
characterised by the 3-Riemann tensor lZ abc d, which when contracted provides the spatial Ricci 
tensor 1Z ab = TZ c acb and Ricci scalar TZ = TZ a a ^\ For our purposes, the key expression is the 
Gauss-Godazzi equation 

TZab = \TZh ab + \li ab - \® (a ab + u ab ) + cr c ( <7 c 6 ) + 2a c [ a uj c b ] + E ab , (26) 

implying that the 3- Ricci tensor has also an antisymmetric part. In fact, lZ ab decomposes into 

TZ{ab) = \TZh ab + \H ab - \®o- ab + <J c ( a a c b ) + E ab (27) 

6 We remind the reader that in rotating spacetimes there are no proper spacelike hypersurfaces orthogonal to 
the fluid flow. Here, we use the term 'spatial', when referring to the instantaneous rest space of a fundamental 
observer, for convenience only. Also, when uo a 7^ 0, H a bcd does not share the full symmetries of the spacetime 
Riemann tensor ES. 
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and 

K a = -\<dui a -a ab J , (28) 

where TZ a = ^e abc lZ bc represents the skew part of lZ ab . The average spatial curvature is de- 
termined by the sign of the 3-Ricci scalar. The latter follows from the generalised Friedmann 
equation 

K = 2(j, (l + \h) - |G 2 + 2a 2 + 2A , (29) 

and evolves according to 



TZ = -§G7£ + 2 



(a 2 ) ' + 2Ga 2 } - 2a ab U ab - ±®A , (30) 



For the Weyl tensor components, we use Eqs. (|24|) and (|2q ) to obtain a set of two propagation 
formulae, namely the E'-equation 

E {ab ) = -@(E ab - in a6 ) + 3 (a c(a + ±e C(% ^) (E c h) - ±II C 6> ) - \^{l + w)a ab 

+(curl#) o6 , (31) 

and the //-equation 

H ab = -@H ab + 3a c{a H c b) - (curl£) a6 + i(curin) afc + 2e cd{a E c b) u d , (32) 
which are accompanied by an equal number of constraints, the (divE')-equation 

V b E ab = J^(v a + \h£ a ) - ±e abc H b curlH c + e abc a b d H cd , (33) 
and the (div/T)-equation 

V b H ab = n(l + w + %h)u a - e abc a b d (E cd + \li cd ) + 3 (E ab - \H ab ) J , (34) 

where B a = SV a H 2 /H 2 . The latter is an additional gauge-invariant variable representing 
inhomogeneities in the magnetic energy density. 

(v) The linear propagation equations of the key inhomogeneity variables V a , Z a and M. a b- 
The former describes gradients in the fluid energy density and evolves according to 

OS 

V a = wQV a - (1 + w)Z a + e abc H b curlH c + § h@Su a 

-a ab V b + U ab ii b , (35) 

A* 

obtained from Eq. (71) in ||. Note the last two terms in the right hand side of the above, which 
reflect the background anisotropy. Following Eq. (72) in ||, the expansion gradients propagate 
as 

Z a = -\QZ a -\iiV a + SA a -3Sa 2 u a + ^Se abc H b cmlH c - ^hB a 

-a ab Z b + ISU ab u b - 2SV a a 2 , (36) 
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where, A a = V a ^4 and the last three terms are due to the anisotropic background. Finally, 
inhomogeneities in the magnetic field vector are governed by 

M ab = ~lQM ab - \K a Z h - \<dS (H [a u b] + 3# (a u b> ) + SH c V b (a ac + e acd to d ) 

-Se acd H c H d b + a ac M% - a bc M a c + SH c u c a ab + 2SH c a c[a u b] , (37) 

which derives from Eq. (75) in ||. Here, the background anisotropy is represented by the last 
four terms of the right hand side. 



4 The Bianchi I Background 

I has long been known that anisotropic exact Bianchi cosmologies of type I can accommodate 



spatially homogeneous ordered magnetic fields. We refer the reader to [20|-|23] for a sample 
of such studies during the 60's and the early 70's. All these authors assume, either directly or 
indirectly, that the magnetic field is aligned along the shear eigenvector, and the same is also 
true for our analysis. Note, however, that the field equations do not necessarily impose any such 
condition upon the magnetic vector in Bianchi I models [24]. In fact a magnetised Bianchi I 



universe, where the field is not a shear eigenvector, has been studied by means of dynamical 
system techniques [p5f|. 

Here, we focus upon the anisotropic magnetic stresses and their impact on cosmological 
perturbations. For our purposes, the whole of the background anisotropy must result from the 
presence of the field. This, in turn, will isolate the effects of the magnetic stresses and facilitate 
their analysis. 

4.1 The Anisotropy 

The key effect of any large-scale magnetic field is the introduction of a preferred direction, which 
is clearly manifested by the eigenvalues of the H ab tensor. According to Eq. (|5|), these eigenvalues 
are — ^H 2 along the magnetic direction and \H 2 in the orthogonal plane. In other words there 
is a negative pressure, a tension, parallel to the field lines. Thus, we see the prominent role 
played by the magnetic lines of force. Each small flux tube is like a rubber band under tension, 
and infinitely elastic [26]. We will return to this property in sections 5 and 7, when discussing 



the field effects on the kinematics of the universe and on the growth of density perturbations. 

Given that (n a ftn a '') 1 / 2 ~ H 2 , the dimensionless energy density ratio h = H 2 /fi offers a 
particularly useful measure of the magnetically induced anisotropy, provided that matter is 
distributed isotropically. Clearly, the degree of distortion inflicted on the background isotropy 
by the magnetic presence depends on the field's strength relative to the fluid, that is on the 
value of h with respect to unity. When the field is relatively weak, that is for /i< 1, one might 
treat the background universe as isotropic to leading order. Strictly speaking however, only 
anisotropic spacetimes can naturally accommodate large-scale cosmological magnetic fields. 

Let us assume that the entire background anisotropy is due to the magnetic presence. This 
scheme may not represent the most general case, but it helps to isolate the effects of the magnetic 
stresses. Therefore, to zero order, both shear and Weyl anisotropics result from the presence of 
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the field. Mathematically, we achieve this by demanding that the background magnetic field is 
aligned in the direction of a shear eigenvector. 

v a bH b = \\H a , (38) 

where V a A = to zero order. Now, we introduce a shear eigenframe with the background field 
vector along its x 1 -axis. In this frame, hereafter referred to as the magnetic frame, both a ab and 
U a b are simultaneously diagonalisable (recall that H a is the natural eigenvector of Ti a b). Axial 
symmetry then implies 

(J ab = -Jptlab , (39) 

thus guaranteeing that shear anisotropies result directly from the magnetic stresses. The above 
resembles the Navier-Stokes equation, so that H 2 /2X may be thought of as the magneto-shear 
viscocity coefficient.^] Following Eq. (|38|), the shear eigenvalues are |A, parallel to the field 
vector, and — ^A orthogonal to it. Also, relations ( |3"8| ) and (^) imply that 

= ^a b v ab = ^A 2 , (40) 

to zero order. 

We monitor the impact of the shear, relative to the average volume expansion, via the 
dimensionless shear-impact parameter £ = A/0. The latter measures the degree of the kinematic 
anisotropy. As an example, consider the expansion rate along the field lines. Then, on using 
@) we find 

e ab n a n b = 1(1 + 206, (41) 

where Q ab = cr a b + |®^a& is the zero-order expansion tensor and n a = H a /VH 2 is the unit 
vector in the direction of the field. Although the shear effect depends on the specific form of £, 
it becomes insignificant when ( < 1. The same condition guarantees that the shear impact is 
also weak normally to the field lines. In fact, result ( |4C| ) ensures that the shear effect, relative 
to the volume expansion, is always negligible provided £ -C 1.Q 

4.2 The Evolution 

In this section we discuss key features of the background Bianchi I spacetime, assuming that 
the magnetic field is the sole source of its anisotropy. Unless stated otherwise, the equations 
presented here are obtained from the linearised formulae of Sec. 3 by keeping their zero-order 
terms only. 

7 The direct dependence of the shear, and subsequently of the electric Weyl tensor, on the magnetic stresses 
holds in the background only. No relation analogous to (|3q ) or (^), between any two ansisotropic sources, is ever 
imposed in the perturbed spacetime. 

8 The overall weakness of the shear is also necessary, if the comoving gradients that describe the inhomogeneities 
are to maintain the desirable geometrical and physical interpretation they have within an almost-FRW universe 
(see Appendix). 
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4.2.1 Energy Conservation 

The fluid and the magnetic energy densities obey the familiar conservation laws. These are 
respectively represented by the equations of continuity 

/i = -(l + w)/i0, (42) 

for the fluid, and 

(H 2 )- = -l(l-0@H 2 , (43) 
for the field. The latter derives from Eq. ( |22"|) on using condition (|3§|). 



4.2.2 Volume Expansion 

To zero-order, result fl40| ) simplifies Raychaudhuri's equation (|l^) to 

9 = -i(l + 2C 2 )e 2 -^(l + 3u; + / l )+A, (44) 

where the magnetic and the shear effects, respectively represented by h and £, assist the matter 
in focusing the fluid flow-lines. Similarly to the FRW models, we may introduce an average scale 
factor (S) via S/S = which also defines the mean Hubble rate. 

4.2.3 Shear Anisotropies 

The spatial flatness of the Bianchi I spacetime and the Gauss-Godazzi equation (see Eq. (83) 
in ||) provide an extra evolution equation for the background shear 

<?ab = -QcFab + ^ab , (45) 

thus supplementing the standard zero-order propagation formula for the shear. The latter is 
obtained from Eq. (|j~4|) by means of condition (|39|) and has the form 



a ab = -l(l + \t)<da ab + \Tl ab -E ab , (46) 

which guarantees that the zero-order shear evolution is shaped by magnetic anisotropies and by 
long-range gravitational forces. The flatness of the unperturbed spatial sections and Eq. ( |29{ ) 
also supply a zero-order expression for the magnitude of the shear tensor 

a 2 = -n(l + \h) + ±G 2 -A, (47) 

which is simply the background Friedmann equation. Note that on using Eq. (|4(]) the above 
writes as 

+ -ie 2 (l-C 2 ) +A = 0, (48) 

an expression that will prove useful later. Finally, employing (|3^) and (^), the zero-order 
component of Eq. (|30| ) gives 

& = -@a-^fxh, (49) 

which, on using (flop, is recast as 

X = -@\-lj l h, (50) 
that monitors the evolution of the A-parameter. 
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4.2.4 Magnetic Stresses 

Equations (^) and ( f46[ ) combine to provide the zero-order expression 

Ii ab = \{l-QQa ah -2E ab , 



(51) 



which is also obtained by substituting (|39|) into the background component of Eq. (|23|), The 
above provides a unique relation between H ab , a ab and E ab . This fine-tuning of the unperturbed 
anisotropic sources results from the spatial flatness of the Bianchi I model and ensures that Eqs. 



D are simultaneously satisfied. Using condition (|39|), Eqs. (p4) gives 

ti ab = -l(i-C)en ab , (52) 

to zero order. The latter is also obtained from Eq. ( |39|) by means of (|39|), (43) and fl49|). 



4.2.5 Weyl anisotropies 

In the unperturbed universe, Eq. (51) guarantees the direct dependence of the Weyl anisotropies 
on the magnetic stresses. In fact, Eqs. (38) and ( |5l"| ) confirm that H a is also an eigenvector of 
E ab , while Eqs. (3£) and ([51]) ensure that 



E, 



ah 



H 2 



n. 



ab ■ 



(53) 



where rj = ^nh + ^(l — ()XQ has the dimensions of E ab . The background evolution of the electric 
Weyl tensor is given by the zero-order part of Eq. (pTl) 



E 



(ab) 



© [Eab - 2^-o.b J + 3<7 c ( a (-E C 6) 



it 



i//(l + w)a ab 



(54) 



Hence, condition ( |38|) and axial symmetry guarantee the direct dependence of the whole back- 
ground anisotropy on the magnetic presence. Once the parameter £ has been established, the 
zero-order evolution of both a ab and E ab is determined by Eq. ((BT 

4.3 The Alfven Speed 

Disturbances in a magnetised medium propagate via the Alfven speed. For an arbitrarily strong 
magnetic field we define 



h 



1 + w + h' 



(55) 



so that c 2 < 1 always. This coincides with the definitions in |g] and |1C] for a weak field in 



a dust-dominated universe (i.e. h <S 1 and w = 0), but deviates when the fluid has non- 
vanishing pressure. The reader must take this difference into account when comparing our 
pre-recombination equations to those in || or [10]. 
On using Eqs. 



and (42), we obtain the evolution law for the ratio between the magnetic 
and the fluid energy densities 



-Ul-3w-AQeh. 



(56) 
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Note that when w 7^ 1/3 and the field is weak (i.e. £ <C 1), the last term in the right hand side 
of (|56| ) is negligible. Hence, in the dust era, h evolves exactly as in weakly magnetised FRW 
cosmologies. In the radiation era, however, w = 1/3 and therefore h = AXh/3. In this case the 
behaviour of h depends on the specific form of A, that is on the relation between a ab and H a b- 
This complication, which was first pointed out in ^7|, corresponds to the 'critical case' discussed 

The Alfven speed also shows the same subtle behaviour during the radiation epoch. Indeed, 
combining definition (|55|) with Eqs. (O) and (|56|) we obtain 



(el} = -| (l - 3c s 2 - 4C) (l - cl) Bel • (57) 

Clearly, when c 2 = 1/3, the evolution of the Alfven speed couples to the shear impact parameter, 
even for c 2 <C 1. 

5 The Case of a Barotropic Fluid 



The set of equations (|35|)-(|37|), which holds for a general perfect fluid with pressure, does not 
contain an evolution formula for the pressure gradients. The reason being that the propagation of 
Y a is determined directly form (^), once the material content of the universe has been specified. 
Here we consider a universe filled with a single barotropic perfect fluid. This has the simple 
equation of state p = p(n) which implies that 

SY a = fic 2 V a . (58) 

5.1 The Scalar Variables 

The gauge- invariant variables defined next, allow to distinguish between inhomogeneities of 
different types. To begin with, we decompose the comoving spatial gradient of T> a as 

SV b V a = A ab = lAh ab + e ahc W c + S afc , (59) 

where A = A a a , W a = ^e abc A bc and T< ab = A/ ab \. The scalar A describes the average gravita- 
tional clumping of the fluid, the vector W a represents rotational instabilities and T, ab describes 
volume-true infinitesimal distortions in clustering |Tj]_ \y e a i so apply similar decompositions 
to Z ab = SV b Z a , for the expansion inhomogeneities, and to B ab = SV b B a , which is associated 
with gradients in the magnetic energy density. 

Given that density perturbations (i.e. scalar matter aggregations) are the focus of this 
analysis the following scalars are crucial 

A = A a a = ^vV Z = Z a a = S 2 V 2 e and B = B a a = ^V 2 H 2 , (60) 

where the equalities hold in the linear regime only. The above variables, which are all gauge- 
invariant, are supplemented by the first-order scalar 

K = S 2 1Z, (61) 

representing spatial curvature perturbations. Note that only A, B and fC are dimensionless, 
while Z has inverse time dimensions. 



14 



5.2 The Kinematic Evolution 



5.2.1 Acceleration 

Employing the barotropic relation (|58[), the momentum density conservation law (see Eq. (|i~0|) ) 
gives 

with the background anisotropy reflected in the last term on the right. Now, expression (||) 
recasts Eq. ( |o2| ) into 

u " = - TTT W i 1 ~ c ") V - ~ 4-Za b cH b curlH c - n i C \ Q n a n b V b , (63) 
(l + w)S \ ' \m [l + w)b 

where the last term acts only along the direction of the field and the second last is transverse. 
Hence, the acceleration components are 



(i+w)s ' Da parallel to H a , 

(64) 

1 - eg) V a - §-e abc H b cmlH c orthogonal to H a , 



(l+w)S \ ^> a iih 

confirming that the field does not affect the fluid motion along its lines of force. Note that, when 
the field is weak (i.e. for eg <C 1) the acceleration components reduce to their FRW counterparts 
(compare to Eq. (31) in |i~0|| ). 

Equation (p2j verifies that, in the absence of fluid pressure (i.e. when eg = 0), ii a depends 
only on magnetic gradients and is always normal to the field vector. Clearly, this is not the case 
when p ^ 0. For non-vanishing pressure, one might also say that H a is no longer normal to the 
fluid flow. We verify these statements by simply contracting (|6^ ) with H a . We find that 

H a u a = H a u a = - -a , c H a V a , (65) 
(1 + w)S 

where generally H a D a ^ 0. The above condition, which is independent of the field's strength, 
is identical to that obtained by the almost-FRW treatment p 



According to the kinematic equations (p^)-(14), the acceleration gradients affect all aspects of 
the expansion. They also have an impact on the spatial curvature (see Eq. (p0|)). Consequently, 
the following linear decomposition is particularly useful 



+H^ab + h / s ^ w2 n a c A cb - ^(1 - C)^a c e bcd u; d - ^-Tl cd K acbd , (66) 
° (j,h{l + w)S z 3/x/i /in 

where the last three terms result from the anisotropic background. Note the magneto-curvature 
contributions in the right hand side of the above. Their presence manifests the vectorial nature 
of the magnetic field and underlines the coupling between magnetism and geometry. Using Eq. 
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(Eft, the two magneto-curvature terms reduce to one given by clH c H d TZ acbd /l^h. The latter 
represents a force transverse to the background field lines, as the symmetries of TZ abc d confirm. 
This term, which appears inevitably when the 3-gradients of the field vector commute (see Eq. 
(B9) in fT6|), can trigger a range of rather unexpected effects. 

Next the trace, the skew part and the symmetric trace- free component of Eq. (66) will be 
used to analyse the magnetohydrodynamical effects upon the kinematics, the dynamics and the 
geometry of our cosmological model. 



5.2.2 Deceleration Parameter 

On using expression (g), the linearised trace of fl6q) gives 

'2 



A 



(l + w)S 2 



A ab n a n b + (l - cfj A ai 



ab 



-^B + clK ab n a n b , 



(67) 



where A = V a u a to first order (see Eq. (|i~2|)) and v ab = h ab — n a n b projects into the 2- 
plane orthogonal to H a . The scalars A ab n a nr and TZ ab n a n b respectively represent density and 
3-curvature perturbations parallel to the background field lines. On the other hand, A ab v ab 
describes matter aggregations orthogonal to H a , so that A ab n a n b + A ab v ab = A. Note that 
S ab n a n b + S ab v ab = S a a for any spacelike tensor S ab , while S ab n a n b = S\\ and S ab v ab = S^ + S^ 
in the magnetic frame. 

Substituting the above result into Eq. (12) we obtain the Raychaudhuri equation for a 
magnetised almost Bianchi I universe filled with a perfectly conducting ideal fluid. The resulting 
formula is then transformed into the following alternative expression 



I6 2 q = k(l + 3u) + A) + 



(1 + w)S 2 
c 2 K ab n a n b + 2a 2 - A , 



A ab n a n b + (l - c 2 ) A ab v ab 



+ 



2S 2 



B 



(68) 



where q = —SS/S 2 is the dimensionless deceleration parameter. Clearly, the sign of the right 
hand side determines the state of the expansion, with positive terms slowing it down and negative 
accelerating it. As expected, the field slows the universe down by adding to the total energy 
density, and the shear has a similar decelerating impact. Comparing to Eq. (36) in |10| |, we 
see that the effect from local increases in the magnetic energy density (i.e. for B > 0) remains 
unchanged. However, the anisotropic background has left its signature through a number of 
direction dependent effects. In particular, for relatively strong magnetic fields (i.e. when c 2 ~ 1), 
matter aggregations contribute differently along the field direction than perpendicularly to H a . 
In fact, as the magnetic field gets stronger (i.e. as c 2 — > 1), the decelerating effect from A ab v ab 
tends to zero. 

According to Eq. (p8|) , only curvature perturbations parallel to the magnetic force- lines have 



an influence on the expansion rate. This effect was originally identified in [10|, although there, 
the isotropy of the background FRW model meant that it was independent of direction. Here, 



as in [10], negative curvature perturbations add to the total deceleration, whereas a positive 
lZ ab n a n b tends to accelerate the universe. This sounds odd at first, given that positive curvature 
is traditionally associated with gravitational collapse. However, the tension carried along the 
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direction of the field, means that small magnetic flux tubes behave like elastic rubber bands 



1 26]. When coupled to geometry, this property of the field gives rice to a curvature stress that 
opposes the action of the magnetic energy density gradients. Hence, in Eq. fl68]), the magneto- 
geometrical term reacts to local changes in the spatial curvature by modifying the expansion rate 
of the perturbed region accordingly. Intuitively, the coupling between magnetism and geometry 
has injected the elastic properties of the field into space itself. Note that the effect of the 
relativistic magneto-curvature term c 2 lZ a bn a n b in Eq. (|68|), bares a striking resemblance to 
the classical (non-relativistic) curvature stress exerted by a distorted field. Following p6[j , the 
tension from field lines with a local radius of curvature R, exerts a transverse force ~ H 2 /R per 
unit volume. Here, the distortion in the field pattern is triggered by perturbations in the spatial 
geometry itself. 

Comparing to Eq. (36) in [10|, we note that the magneto-curvature effect identified in Eq. 



reduces to that predicted by the FRW-based analysis when TZ a i,n a n b = \TZ. Thus, for a 
weak field, the two equations become indistinguishable when the curvature perturbation along 
H a takes the (directional) average value. Generally, however, lZ a bn a n b has a rather complicated 
evolution that also depends on magnetic, shear and Weyl contributions. In fact, contracting Eq. 
( |26| ) with the field vector we find 

n ah n a n h = \K - Ifih - fcr 2 - |(1 - 2C,)Q<j ab n a n b + E ab n a n b , (69) 

where a a \,n a n b and E a t,n a n b are respectively the components of the shear and the electric Weyl 
tensors parallel to the filed lines. 

The anisotropy of the Bianchi I background has introduced directional effects to the expan- 
sion dynamics, along and normally to the background magnetic field lines. As we shall show next, 
similar direction-dependent complications affect almost every aspect of the evolution. None of 
these effects could have been identified by an FRW-based treatment, where the isotropy of the 
unperturbed spacetime allowed only the average of these effects to be measured. 

5.2.3 Vorticity 

Substituting the background relation ( |39| ) to the vorticity propagation equation (|i~3| ) we have 

w« = -§ (l + |c) ©^a - |curlu a , (70) 

where £ always represents the magnetically induced kinematic anisotropies. Using Eqs. (||), 
(|28|), (|38|), (|3£), (f43|) and (|59|), the commutator (B8) from |l6|| , the symmetries of 1Z a bcd (see 
eq. (81) in p8|) and the linear relation W a = —(1 + w)@S 2 ui a , we obtain an expression for the 
skew part of Eq. (^). Then, given that curlu a = e a bc^ h u c , the vorticity propagation equation 
is recast as 



(l + K) " 1 (l " VS) i] = -^~ h H'WH a+ 2lih{ ff m)s ^ ^ d , (71) 



where X a fc = (S 2 / (i)V ( a V ^ to first order. The above holds for an electrically neutral medium 
with H a uj a = 0. As predicted by the almost-FRW treatment (see Eq. (37) in [|10|), magnetism 
becomes a source of vorticity when cvLi\H a varies in the direction of the field vector. Note 
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that before recombination distortions in the density distribution can also generate vorticity. 
According to Eq. (fn|), a relatively strong magnetic field with £, ~ 1 can modify the vortex 
dilution caused by the volume expansion. 

5.2.4 Shear 

The symmetric trace- free part of (66) transforms Eq. (|14|) into 

c 2 c 2 

r 2 ~ f? 



Comparing with Eq. (38) in p0| , we identify the effects of the background anisotropy in the last 
four terms. The extra magnetic input to the S a fe-term, which is negligible for relatively weak 
fields, results from the magnetic contribution to the total energy density of the universe. Any 
other changes are due to differences in the Alfven speed definitions. 

An additional useful relation is the propagation equation of the shear magnitude. We obtain 
it by contracting Eq. (83) in |Bj with a a b- We then employ results ( |66[ ) and (|67|), the symmetries 
of IZabai and the constraint M. a a = to find 

a 2 ) ' = -20a 2 - ^ha ab n a n b - ^C© [2 (l + \cl) K ab n a n b - K ab v ah 



S<;B [ 2A afe nV - (l - cl) A ab v' !h 



3(1 + w)S 2 



{; s . 2 (2B ab n a n b - B ab v ab ) , (73) 
where a ab n a n b satisfies Eq. (|69|). Also, B ab n a n b and B ab v ab respectively represent magnetic 



energy density gradients parallel and orthogonal to the background field lines. Similarly, TZ ab v ab 
is the 3-curvature perturbation in the plane normal to H a . Note that the anisotropic effects, 
manifested by the direction dependent terms of the right hand side, remain even when the field 
is relatively weak. The above result will prove useful in Sec. 5, when studying spatial curvature 
perturbations. 

5.3 The Dynamic Evolution 

5.3.1 Perturbations in the Matter Density 

In the linear regime, magnetised density gradients propagate according to Eq. (35). For a 
barotropic perfectly conducting medium, the latter becomes 

Va = (w- c 2 cl + ±C) ®V a - (1 + w)Z a + ^ (l - cl) e abc H b cnv\H c 

-(c 2 s cl-C)en a n b V b , (74) 
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on using (|5|), the zero-order relation (|3^), the barotropic expression ( |6^ ) and result (|65|). We 
can now see that the field affects density perturbations along as well as normally to its own 
force- lines. Parallel to H a , the density gradients are influenced by the magnetically induced 
kinematical anisotropics, with their presence manifested by the shear impact parameter £. 

The evolution of linear matter aggregations is determined by the 3-divergence of Eq. (|74|). 
Using the commutation law between time derivatives and 3-gradients (see formula (B18) in |16|]), 
we obtain 

A = (w — f C) &A ab n a n b + (w - c 2 c 2 + |c) &A ab v ab - (1 + w)Z 

+ + w)cleB - (1 + w)cleJC ab n a n b . (75) 

Note that different quantities contribute differently parallel than perpendicular to the field lines. 
Also, the magneto-shear effects, which are imprinted in the impact parameter vanish when 
A ab n a n b takes the directional average value |A. 



5.3.2 Perturbations in the Expansion 

Starting from Eq. (|36|) we employ expression (||), the zero-order relations (^), (|38| ) and (|47|), 
together with the linear equations (|62[), (|65|) and (37) to arrive at the propagation formula of 
the expansion gradients in a barotropic fluid environment 



(i - |c) ez a - \ 



l + 4c 2 C 2 (l + §r) 



lnhB a + | 



1-fc: 



a (l + §r) 



Se abc H b cmlH c 



(l + w)S 



A ab n a n h + (l - cl) A a 



ab 



bV 



c ~ 
— V B 



<Qn a n b Z b + ^ h ° l 



1-W 



l + |r 



2(1 + u> 

+c 2 SV a (ft ab n a n b ) - 2C95V a (a ab n a 



n a n V b 



(76) 



where the dimensionless parameter r = ^ 1 — ^ ^|0 2 — A^) will prove useful later. 

Taking the comoving spatial divergence of Eq. (f76|), using the commutator between time 
derivatives and spatial gradients of first order spacelike vectors (see Eq. (B18) in [|l(|), together 
with expression (0) and the background relation (pjJ) we find that, to first order, 



Z 



2-e 



(1 + 2QZ ab n a n b + (1 - ()Z ab v 



ab 



1 + 



1+W 



(1 +6r) 



a„b 



A ab n a n 



l + 4c 2 c 2 (l + |r) 



06 



1 - 4c 2 (l + |r 



B-\iih 



l + |r 



K, ab n a n b 



t , . A ab n a n b + (l - c 2 ] A a , 
+c 2 V 2 (/C^nV) - 2(QS 2 V 2 (a ab n a 



ab 



¥lV 2 B 



(77) 
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where Z ab n a n and Z ab v are respectively the expansion gradients parallel and orthogonal to 
the zero-order magnetic vector. Note that most of the anisotropic effects in (|76|), (f77| ) die away at 
the weak magnetic field limit. Those that remain are the directional dependence of the curvature 
contribution, together with the curvature and shear Laplacians. 



5.3.3 Perturbations in the Magnetic Energy Density 



For the linear evolution of the magnetic energy density gradients we start from the definition of 
B. Taking the time derivative of Eq. (|60|c)), we employ expression (||), the linear conservation 

(guaranteing that a ab n a n b 



law (22), the background relations (| 
the linear result (B7|) and the propagation equation (fT5|) to find 



| A to zero-order) and 



B 



3(1 + w) 
46 



3(l + u;) 



4@ 



|[ W _ C 2 ( 1_ C) _4 C 



C (i + U 2 



2^a 



B ab n a n b - |C (l - cl) B a , 



A ab n a n b + 



w 

ab 



l-C(l-cl)) +§C] A ab v ab } 



+jcl(QlC ab n a n b + 2S 2 V 2 (a ab n a n l 



(78) 



where the shear Laplacian is the only effect of the anisotropic background that remains at the 
weak magnetic field limit. In deriving Eq. (|78|) we have employed the linear commutation law 



(vv)- 



v 2 / 



\eV 2 f-2a ab V a V b f + fA, 



(79) 



where V a / vanishes in the background. The above follows from linearising the 3-divergence of 
(B18) in |1§]. 



5.3.4 Perturbations in the Spatial Curvature 

On using the barotropic expressions (67) and (f7|), Eq. ( ^0|) is recast as 

U = -|G{[l + 2c 2 (l + ic) +2c\n ab n a n b + {l-QK ab v ab ) 



+ 



4c 2 e 



+■ 



3(1 + w)S 2 

2de 



(l-C)A afe nV+(l-c 2 ) (l + K) A < 



ab 



35 2 



(l-()B ab n a n b + (1 + k )BabV 



ab 



(80) 



which reduces to Eq. (40) in [|l(]] at the weak field limit. For a relatively strong field with 
c 2 , ( ~ 1, the kinematic anisotropies (represented by £) complicate the overall magnetic im- 
pact on curvature. One faces a highly anisotropic situation, where disturbances parallel to the 
background field increase the average curvature, while perturbations orthogonal to H a have 
the opposite effect and vice versa. If we ignore all (^-contributions in Eq. (|80|), the magneto- 
curvature effect on 1Z is imprinted in the direction dependent term 2c 2 JZ ab n a n b . The latter 
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boosts the smoothing effect of the expansion on 3-curvature perturbations in agreement with 
[10]. In fact, the two effects coincide when lZ a i,n a n b = ^TZ. 

It should be emphasised that all our equations refer to an arbitrarily strong magnetic field. 
So far, the only restriction is that the background field is an eigenvector of the shear tensor, 
which translates into a simple algebraic dependence of the zero-order shear and electric Weyl 
tensors on li^ (see Eqs. ( p9| ) and (|53|)). Next, we will apply these formulae to the case of a 
relatively weak background magnetic field. 



6 The Weak Magnetic Field Limit 

The observed anisotropy of the CMB imposes strict limits on the strength of any potential 
primordial magnetic field at the time of decoupling. Also, Helium-4 measurements restrict the 
magnitude of such field at the time of nucleosynthesis. All point towards a relatively weak 
magnetic field, with current estimates arguing for a field strength of ~ 10~ 8 G at present. Here, 
on the basis of these observations, we assume a background magnetic field weak compared to 
the dominant isotropically distributed matter component (i.e. h, c 2 ~ h/(l + w) <C 1). Given 
that magnetism is the sole source of the background anisotropy, it is reasonable to assume that 
the shear impact is also weak to zero order (i.e. £ <C 1). At this limit, the unperturbed and 
linear equations simplify considerably. 



6.1 Background Evolution 

Given that the source of the anisotropy is weak, one expects that certain key figures of our 
Bianchi-I background will approach their Friedmannian counterparts. Indeed, when £ <C 1 the 
magnetic energy density (see Eq. fl43|)) evolves unaffected by the background anisotropy 

ff 2 )' = -§0# 2 . (81) 

The same is also true for the fluid energy density, given the absence of magnetic terms in Eq. 



The weakness of the field also means that the average volume expansion of our Bianchi I back- 
ground approaches that of the Friedmann model. When h, £ <C 1, the zero-order Raychaudhuri 



equation ([44|) gives 

9 = -±6 2 - i/x(l + 3u>) + A, (83) 

as in isotropic FRW cosmologies. The generalised Friedmann equation ( |48| ) also acquires an 
FRW profile 

//-±G 2 + A = 0, (84) 

9 For radiation w = | and Eq. (^6|) implies that h — j-Q£h. Also, the shear impact parameter evolves as 
£ = — |0(C + h), given that £ = A/0 and using Eqs. (|50|), (pS|), ( p4| ) and A = 0. Thus, when there is no magnetic 
field (i.e. h = 0), the expansion anisotropy drops rapidly. The field preserves the anisotropy so that £ does not 
tend to zero. Instead, £ — *• at late times and h experiences a slow logarithmic decay 

k = l + 2Mi/t ) ' (82) 
in agreement with the 'quasi-static' solution obtained in MM. 
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given the spatial flatness of the Bianchi I spacetime. This guarantees that the dimensionless 
parameter r in Eqs. (|76|), ( |77| ) vanishes at the weak field limit. 

It should be emphasised that these results do not, in any case, suggest that the unperturbed 
Bianchi I model has now been reduced to a simple FRW cosmology. Our background universe 
is still anisotropic, supplied with magnetic stresses, shear and an electric Weyl tensor. However, 
the weakness of the magnetic field means that, to leading order, the anisotropy no longer affects 
certain key evolutionary aspects, such as the rate of the average volume expansion. 



6.2 Linear Evolution 
6.2.1 Kinematics 

For a weak background magnetic field, Raychaudhuri's equation, written in terms of the decel- 
eration parameter, becomes 

2 2 

±G 2 q = \n{l + 3w) + 2a 2 + ^' A + - c\ll ah n a n h - A . (85) 

(1 + w)b z lb z 

Relative to Eq. (|68|), all the directional dependences in the contributions from matter aggre- 
gations have been smoothed out. On the other hand, comparing to Eq. (36) in [ 10 1 , we notice 
the extra shear term, which adds to the total deceleration. The anisotropic background has also 
refined the curvature contribution, which is now confined along the field lines only. In agreement 
with the FRW treatment, the magneto-geometrical term accelerates models with positive local 
curvature, but slows down negatively curved regions. The effect observed here reduces exactly 
to that identified by the FRW analysis when lZ a bn a n b = ^1Z. Generally, however, 

1l ab n a n b = \K- l(J,h - fcr 2 - \@cr ab n a n b + E ab n a n b . (86) 

Note that the impact of the <r 2 -term in Eq. (^5|) becomes negligible provided that the weak field 
restriction is imposed upon the actual universe, as well as the background. 

Given the weakness of the background magnetic field, the vorticity propagation equation 
reduces to 

*. + §(!- K) e u „ = -^l—^H'-v^ + ^ff^ ^n^ . (8 7) 

For dust, the above agrees with the FRW equation (37) in @. 

Assuming a weak background magnetic field has very little effect on the shear evolution. In 
fact, for c 2 , C < 1 Eq. fl7|) gives 

c 2 c 2 



T(ab) 



+ ^hS H ^ cM ^ - Eab - + ^h(l S +w)S^ aACb 

r 2 dr 2 f) 
a tt t> c d *a" n c ^ , A 

-— u cd n {ab) - —n {a e b)cd u , 



which is considerably more complicated than Eq. (38) in [10|. This is not surprising, given that 
the shear reflects the intrinsic spatial anisotropy of the Bianchi I spacetime. 



22 



6.2.2 Dynamics 

At the weak field limit, Eq. (f7^) that monitors gradients in the fluid energy density, reduces 

ton 

A = wOA - (1 + w)Z + |(1 + w)clQB - (1 + w)cl®K ab n a n h , (89) 
Again, all directional dependences in the contribution of A have died away. The difference 



between the above and its FRW counterpart, given by Eq. (54) in [10], is in their magneto- 
curvature terms. The two equations coincide when JC ab n a n b = |/C. 

When the field is relatively weak, r = and Eq. (f77|), for the propagation of the expansion 
gradients simplifies to 

Z = -|GZ-|/xA + i(l + u;)/xc^-|(l + y;)/ic^/C Qb nV 

2 

V 2 A - \clV 2 B + c 2 V 2 ()C ab n a n b ) - 2(OS 2 V 2 (a ab n a n b ) , (90) 



since all the directional differences in the Z, T> and B contributions have vanished. Nevertheless, 
the anisotropy of the Bianchi I background is still imprinted in Eq. ( |90[ ) via K, ab n a n b and the 
new shear term, both of which are direction dependent. Note that the Laplacian contributions 
become negligible on wavelengths larger than the Hubble horizon. Also, when JC ab n a n b takes 
the average valuei/C all the differences in the spatial curvature terms, between Eq. (^) above 
and Eq. (55) in |L0|], disappear. 

At the weak field limit, almost all the directional effects in Eq ([78]) are removed. The 
magnetic energy gradients propagate as 

3(1 + w) 3(1 + w) v / 

where terms quadratic in the smallness parameters c 2 and Q have been left out. Relative to Eq. 
(61) in [p.C[|, there is an extra input depending on the shear component parallel to H a . This 
effect, which results from the background anisotropy, becomes negligible on superhorizon scales. 

Finally, a weak background magnetic field means that the all anisotropic effects in Eq. 
are smoothed out. The latter reduces to its FRW counterpart 



4c s 2 e 

3(1 + ™) 



1C= „ s , ^ A + ^cleB, (92) 



where 3-curvature perturbations are now represented by the dimensionless scalar /C = S 2 1Z 
(compare to Eq. (57) in p 



10 Looking at Eq. (fTq), it becomes clear that the (^-contributions along and transverse to the field lines are 
relatively weak when w, 7^ 0. Even in the dust era, when w = = c^, the effect of these terms is negligible. 
Indeed, on using definition ( |5S| ) one can verify that A ± ^"QA ~ A for £ <C 1. Similarly, one shows that 
A ± Ca9A ~ A when cl < 1. 
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7 Particular Solutions 



7.1 Large scales 

On wavelengths larger that the Hubble horizon the Laplacian terms in Eqs. ( p0|) and (^Tj) are 
negligible. On these scales, the evolution of magnetised density gradients is monitored by the 
following system of ordinary differential equations 

A = w6A- (l + w)Z + \{l + w)cl®B - {l + w)cl@lC ab n a n h , (93) 

Z = -§OZ -^A + \(l + w)ficlB-^(l + w)ficl)C ab n a n b , (94) 
4 ■ 4 (c 2 - w) © 



* = «frb A + Hes. (96) 



When JC a i ) n a n b = the above coincide with the equations obtained by the FRW-based analysis 
of |l(J.[[]] In this occasion, which one might call the case of 'average' or 'isotropic' curvature 
perturbations, the predictions of the Bianchi I analysis are identical to those of the FRW-based 
treatment. 

Given that the above system closes whenever JC a i,n a n b oc /C, it is worth considering two addi- 
tional critical cases. The 'minimum' curvature case, when JC a bn a n b = 0, and that of 'maximum' 
curvature with JC a bn a n b = JC (i.e. for JC a bV ab = 0). In the former situation, the equations are 
clear of any magneto-curvature complexities and the evolution of A proceeds as predicted by 
the FRW treatment. In the case of maximum curvature, however, there is a small numerical 
difference between the Bianchi I equations and those of the almost-FRW analysis. Note that 



the isotropic and the zero curvature cases have been already addressed in [1C] and subsequently 



refined in [11]. Here, we will reproduce these results for comparison reasons and concentrate on 



the case of maximum curvature. Despite the fact that all three cases are rather special, they are 
quite important. As we shall see next, their study illustrates the effect of the magneto-curvature 
coupling on the growth of linear density perturbations. 

7.1.1 Radiation Era 

During the radiation era w = c 2 = |, the fluid density falls as \i oc S~ A and the Alfven speed can 
be treated as a constant (i.e. (c 2 )' = 0).f^\ Expressed in terms of the scale factor, to facilitate 
the solution, the system (|9^)-(p6D writes 



S 2 ^^ = 2AM _ \ C \B^ + 8#C("> , (97) 
do 



n Any numerical differences between our equations and solutions and those of |10L [fj] are due to differences in 
the Alfven speed definitions. 

12 When radiation dominates the evolution of the background Alfven speed depends on f, the shear impact 
parameter (see Eq. (EtI)). At the weak field limit this dependence translates into a non- linear logarithmic decay 



for cj, identical to that of h (see Eq. (|82])). In any case, the coupling between c 2 and /i or in Eqs. Q93[) , (| 
and (|96h, ensures that the Alfven speed behaves as if it was time independent. 



24 



S 



/CM 
"dfT 



A^ + 2c?flM . 



(98) 



(99) 



where the expansion gradients have been decoupled from the equations and we have set A = 0. 
Also, A( v ) is the z/-th harmonic component of A and similar decompositions have been applied 
to the rest of the inhomogeneity scalars. The above has the following power law solution 



(100) 



where Aq , A^ ar e arbitrary positive constants and the parameter a satisfies the cubic equation 



a 



a 



1 - ^c 2 



a 



8ct 



(101) 



to lowest order in c 2 . This has one positive and two negative roots that correspond to one 
growing and two decaying modes. For small c 2 , the roots are perturbatively found to be 



a 



-l + fcl + 0(ci), 
2 + ld + 0(4) , 

-4*2 + 0(4), 



(102) 



which implies the following evolution 



AM = + Af } 5 2+ 9 c " + A^S- 1 



-f^ + At^"^ 



(103) 



for the density contrast. In the absence of the magnetic field, namely for c 2 = 0, the above 
reduces to the standard solution of cosmological density perturbations (see [32] for example). 
Hence, the field presence has led to a slight increase of the growing mode. As we shall see 
next, this increase results from the coupling between magnetism and geometry. In particular, 
for minimum curvature perturbations, JC a i,n a n b = and 



a 



-l + lcl + 0(ci) 



(104) 



to lowest order in c 2 ficfl . Thus, in the absence of any curvature contributions the pure magnetic 
effect is to reduce the growth rate of the density contrast. After recombination, an analogous 
reduction of the growing mode was also identified both by the relativistic analysis (see |lOt| ) 
and by the Newtonian treatment (see Q). In [O], this particular effect was attributed to the 
'frozen in' property of the field. The latter combines with the dilution of the magnetic force-lines 
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caused by the expansion to oppose gravitational collapse. Alternatively, for isotropic curvature 
perturbations K. a bn a n b = and 

-1 + 1^ + 0(4), 

a - { 2 + 0(4), (105) 

to lowest order in c\ [|ll]] . In this occasion the density contrast grows exactly as in non- magnetised 
cosmologies. 

Comparing results dT02|) , (|o|) and ( gOp we notice that as the magneto-curvature contribu- 
tion gets stronger, from JC ab n a n b = to K, ab n a n b = |/C and finally to JC ab n a n b = )C, the growing 
mode of A successively increases. In fact, by the time the curvature input has taken its maximum 
value, the pure magnetic effect has been completely reversed. In this case, the density contrast 
grows faster than its counterpart in non-magnetised cosmologies. Again, the reason behind this 
effect is the special form of the anisotropic magnetic stress tensor. The latter ensures that there 
is a negative pressure, a tension, along the direction the field lines. As mentioned in Sec. 5, 
when combined with geometry, this property of the field opposes the action of the magnetic 
energy density gradients (see Eqs. (|93|) , (|97|)). 

Undoubtedly, the weakness of the background field means that the aforementioned magnetic 
effects are small 'corrections' to the standard evolution of density perturbations. Interestingly, 
however, whether the field boosts or inhibits gravitational collapse may entirely depend on the 
coupling between magnetism and geometry. 



7.2 Dust Era and Small Scales 

In the dust era Eqs. (|93|), (|94|) and (p6[) guarantee that, to lowest order in c 2 , 3-curvature has no 
effect on the propagation of A. In other words, after recombination long- wavelength magnetised 
density inhomogeneities proceed exactly as predicted by the Friedmannian treatment |l^, pj]l . 
Thus, it becomes clear that the FRW framework offers a quite accurate treatment of large-scale 
magnetised cosmological perturbations in almost all situations. 

For a general K, ab n a n b , however, equations (|89|)-(|92|) are not adequate to describe the evo- 
lution of A. In order to close the system, one requires the propagation equation for JC ab n a n b 
instead of (p^). This, in turn, is immediately coupled to the shear and the electric Weyl tensors 
via the first order relation 

IC ab n a n b = \K- \S 2 nh - §SV - \QS 2 a ab n a n b + S 2 E ab n a n b , (106) 

obtained form Eq. (|S6"|). It follows that when the anisotropy in the curvature perturbation is 
'non-critical' the problem becomes very complicated to address analytically. In this well 
as on subhorizon scales, one has to resort to numerical methods for solutions. Alternatively, 
one might employ dynamical system techniques (see |3l|] for an overview) to obtain qualitative 
answers. This approach has been employed by (2[| to analyse an exact Bianci I cosmological 
model permeated by source-free large-scale magnetic field. Dynamical system methods were also 
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used in |l6| to study qualitatively non-magnetised density inhomogeneities within perturbed 
Bianchi I cosmologies. Note that one may invoke fluid backreaction arguments and ignore, to 
first order, the small scale shear effects. In this case the system (p9|)-(p2[) reduces to its FRW 
counterpart on every scale provided K, a bn a n b = ~/C juj. 



8 Conclusion 

We have pursued the study of magnetised cosmological perturbations within a perturbed Bianchi I 
spacetime, in the first attempt of this nature. So far, the subject has been only addressed within 
the limits perturbed FRW cosmologies, given the observed high isotropy of the CMB spectrum. 
Strictly mathematically speaking, however, a spatially isotropic spacetime cannot naturally ac- 
commodate large-scale magnetic fields. Even when the background field is is treated as random 
or weak, there are problems related to the generic anisotropy of the magnetic stresses. Although 
it is physically plausible that weak fields can be adequately treated within the almost-FRW limit, 
mathematically speaking the Friedmann spacetimes are unsuitable hosts for cosmic magnetism. 
In this respect, the spatially anisotropic Bianchi I background provides a better framework for 
studying the magnetic effects on cosmological perturbations. Moreover, such study can be easily 
applied to primordial magnetic fields of arbitrary strengths, if required. 

The search for mathematical consistency was not the only motivation behind this work. The 



relativistic, but almost-FRW, treatments of [1C, 11] suggested a number of unexpected magnetic 



effects on the kinematics and the dynamics of the universe. All of them seemed to derive from 
an intriguing coupling between magnetism and geometry. The potential importance of these 
magneto-curvature effects meant that a new study, this time in a more natural environment, 
was necessary. The Bianchi I model, which has long been known to be the simplest anisotropic 
spacetime able to accommodate cosmological magnetic fields, provided just this environment. 

We began by assuming that the entire background anisotropy was due to the magnetic 
presence, in order to isolate the anisotropic effects of the field. Thus, in the background kinematic 
anisotropies as well as anisotropics in the spacetime geometry are both the result of the magnetic 
stresses. In mathematical terms this is achieved by demanding that, to zero order, the magnetic 
vector is an eigenvector of the shear tensor. Although this scheme may not represent the most 
general case, it is sufficient for studying the lowest order effects of the magnetic stresses at the 
onset of structure formation. Any zero-order anisotropies that are not the result of the magnetic 
presence are irrelevant to our study and therefore ignored. At the linear level, however, the field 
is no longer assumed to be a shear eigenvector. Thus, to first order the shear does not entirely 
depend on the magnetic stresses. Our first step was to linearise the non-linear formulae of 
[]|] about a Bianchi I spacetime, and also provide a couple of new but necessary first-order 
relations. As in |llj, we applied our equations to the case of a barotropic medium and then 
examined in detail the magnetohydrodynamical effects on the kinematics and the dynamics of 
the universe. It was at this level that we were first able to identify the implications of the 
anisotropic background. 

In particular, it became clear that some of the effects of the FRW-based treatment were 
in fact direction dependent. To be precise, a number of quantities were found to contribute 
differently parallel to the field lines than normally to them. For example, for a relatively strong 
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magnetic field, the decelerating effect from matter aggregations was no longer isotropic. We 
found that only density perturbations parallel to the direction of the zero-order field vector can 
actually affect the average volume expansion. In contrast, density gradients perpendicular to the 
magnetic force-lines get weaker as the strength of the field increases. Nevertheless, the isotropy 
of this effect was immediately reinstated at the weak field limit to match the predictions of [pX[| . 

The role of the background field stresses was further illustrated through the magneto- 
curvature contribution to the deceleration parameter. In |To|| , the isotropy of the unperturbed 
FRW spacetime meant that the curvature effects were direction independent. Here, however, we 
found that their contribution is entirely confined along the zero-order field vector. As in [ [Tcfl , 
the resulting magneto-curvature term adds to the deceleration of disturbances with local nega- 
tive spatial curvature, but accelerates positively curved regions. This seemed odd at first, given 
that positive curvature is traditionally associated with gravitational collapse. The explanation 
came from the directional dependence of the curvature effect. Recall that there is a negative 
pressure, a tension, in the direction of the magnetic lines of force. This is another way of saying 
that the field lines are elastic, opposing any attempt to distort their equilibrium pattern. When 
coupled to geometry, this magnetic property ensures that the field will react to any changes in 
the spatial curvature by modifying the local expansion rate accordingly. This relativistic effect is 
closely analogous to the classical curvature stress exerted by field lines with a nonzero curvature 
radius. Note that, qualitatively speaking, the aforementioned effect does not depend on the 
field's relative strength. 

Similarly, the role of the field as a source of universal rotation is also independent of its 
strength. As in [10], we find that when the curl of the magnetic vector varies along the direction 
of the background filed, cosmic magnetism generates vortices in the fluid flow. However, the 
magnetic effect on the evolution of these vortices differs according to the strength of the field. In 
particular, the presence of strong magnetic fields can modify the dilution rate of such rotational 
disturbances. 

Our study also facilitated a direct test for the FRW-based treatments of weakly magnetised 
density perturbations. We can see now the number and the nature of the corrections due to 
the background anisotropy, and identify the domains they mainly affect. Therefore, we are in 
a better position to judge the accuracy of the Friedmannian approximation. For an arbitrarily 
strong magnetic field, we find several directional contributions in the propagation equations of 
the basic inhomogeneity variables. As before, the key direction is determined by the background 
magnetic field. At this level, there is a considerable gap separating the Bianchi I equations from 
those of the FRW approximation. At the weak field limit, however, these differences are re- 
duced dramatically. In particular, we find that only curvature and shear effects propagate along 
the preferred direction of the background field vector. Moreover, after recombination, density 
perturbations proceed essentially unaffected by curvature complications. Also, the shear contri- 
butions, which result from the kinematical implications of the background field, are confined to 
sub-horizon scales only. As a result, the large-scale equations of the FRW analysis are identical 
to those of the Bianchi I treatment when TZ a bn a n b , the curvature perturbation along the field 
lines, takes the 'isotropic' value \lZ. Also, the two approaches are essentially indistinguishable 
in two additional critical cases, when TZ a i,n a n b = and for TZ a ^n a n b = 1Z. In the former case 
the curvature perturbation parallel to H a is said to be 'minimum' and in the latter 'maximum'. 
Conclusively, we feel justified to argue that, at least on super-horizon scales, density perturba- 
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tions are accurately treated by the FRW approximation. Even on small scales, one may be able 
to invoke backreaction arguments and ignore the aforementioned shear effects from the linear 
perturbation equations. 

With regard to density perturbations, the Bianchi I analysis has lead to an additional, this 
time less trivial, result. In particular, let us recall the results of the isotropic and the minimum 



curvature cases mentioned above, which were addressed in J10| , 11]. There, it was found that in 
the absence of any curvature input (i.e. for lZ ab n a n b = 0) the magnetic presence suppresses the 
growth rate of matter aggregations. In contrast, when isotropic curvature contributions were 
allowed (i.e. lZ ab n a n b = |7£), the overall magnetic effect was reduced and the density gradients 
grew as fast as in non-magnetised universes. This pattern was also confirmed here. Moreover, 
for maximum curvature contributions, that is when 7Z ab n°"n b increases from ^7Z to 7Z, the density 
contrast grows even faster than in magnetic-free cosmologies. Thus, the coupling of the field to 
the geometry can actually reverse the original magnetic impact on gravitational collapse. The 
reason behind this effect is again the tension properties of the magnetic force-lines, which leads 
to a curvature stress proportional to c\lZ ab n a n b . The latter tends to counterbalance the action 
of the magnetic pressure gradients, thus reversing the pure magnetic impact on linear density 
condensations. As a result, in the magnetic presence, the density contrast may even grow faster 
than in non-magnetised cosmologies. 
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9 Appendix 

9.1 Interpreting T> a 

Following the connecting vector 5x a between two neighbouring worldlines (i.e. fundamental 
observers) satisfies the restriction 

u b V b 5x a = Sx b V b u a . (107) 

Also, the relative position vector, which connects the same two points P, P' on these worldlines 
at all times, is 

<fa(o) = h ab 5x b , (108) 

and propagates as 

(Sx {a) )' = (a ab + e abc io c + \®h ab ) 6x® . (109) 

Suppose now that /J,(x c ) and [i{x c + <5x( c }) is the fluid density at P and P' respectively, 
relative to a frame at P. Then, to first order, 

/i(x c + 5x {c) ) = fj,(x c ) + (V a /i) 5x {a) 

6=— = X a 5x {a): (110) 
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where 5(jl = fi(x c + Sx^) — fi(x c ) and X a = V a /i//i. Equation ( |109| ) ensures that within a FRW 
spacetime 

fe(a> = S (fe(o>) , (HI) 

where \ Sx^ a ^j is a constant. Therefore, result (|11C| ) becomes 

*=(fo< a >) o X) aj (112) 

where P a = SX a . Thus, in an almost-FRW universe the comoving fractional density gradient 
T> a describes the spatial density variations between two neighbouring fundamental observers. 
However, in a Bianchi I spacetime the background shear is non-zero and 8xu\ no longer evolves 



according to (111). Consequently, in a perturbed Bianchi I universe T> a no longer describes the 
spatial density changes seen by a pair of neighbouring fundamental observers. Obviously, in the 
case of weak background shear c/O <C 1 and the aforementioned complication is negligible, and 
D a regains the desired interpretation. 
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